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Abstract — We propose a novel distributed resource allocation 
scheme for the up-link of a cellular multi-carrier system based on 
the message passing (MP) algorithm. In the proposed approach 
each transmitter iteratively sends and receives information mes- 
sages to/from the base station with the goal of achieving an 
optimal resource allocation strategy. The exchanged messages 
are the solution of small distributed allocation problems. To 
reduce the computational load, the MP problems at the terminals 
follow a dynamic programming formulation. The advantage of 
the proposed scheme is that it distributes the computational 
effort among all the transmitters in the cell and it does not 
require the presence of a central controller that takes all the 
decisions. Numerical results show that the proposed approach 
is an excellent solution to the resource allocation problem for 
cellular multi-carrier systems. 

I. Introduction 

Orthogonal Frequency Division (OFDM) modulation is one 
of the candidate technologies for future generation broadband 
wireless networks. Provided that the system parameters are 
accurately dimensioned, OFDM transmissions are not affected 
by intersymbol interference (ISI) even in highly dispersive 
channels. Moreover, OFDM can effectively exploit the channel 
frequency diversity [1], [2] by dynamically adapting power and 
modulation format on all subcarriers. Orthogonal frequency 
multiple access (OFDMA) is the multiple access scheme 
based on OFDM: each user is allocated a different subset 
of orthogonal subcarriers. When the transmitter possesses full 
knowledge of channel state information, the subcarriers can be 
allocated according certain optimality criterion to increase the 
overall spectral efficiency, exploiting the so-called multiuser 
diversity. Resource allocation is one of the most efficient 
techniques to increase the performance of multicarrier systems. 
In fact, propagation channels are independent for each user and 
thus the sub-carriers that are in a deep fade for one user may 
be good ones for another. Many resource allocation algorithms 
have been designed for taking advantage of both the frequency 
selective nature of the channel and the multi-user diversity. In 
most cases dynamic resource allocation has been formulated 
with the goal of either minimizing the transmitted power with 
a rate constraint [3], [4] or maximizing the overall rate with 
a power constraint [5], [6]. 

In this paper, starting from the formulation of resource 
allocation problem as a minimization problem, we propose a 
novel distributed resource allocation scheme for the up-link of 
a cellular multi-carrier system based on the message passing 
(MP) algorithm. MP algorithms have gained their momentum 



in the last years owing to their broad usage in LDPC and turbo 
channel decoding applications [7]. In this setting, messages 
represent probabilities or beliefs 1 , which are exchanged with 
the goal of achieving an optimal bit decisions. We will show 
that resource allocation may rely on a similar MP procedure: 
with the goal of achieving a global optimal assignment, each 
transmitter iteratively sends and receives information messages 
to/from the base station until an allocation decision is taken. 
The exchanged messages are the solution of small distributed 
allocation problems. To reduce the computational load, the 
MP problems at the terminals follow a dynamic programming 
formulation. The advantage of the proposed scheme is that it 
distributes the computational effort among all the transmitters 
in the cell and it does not require the presence of a central 
controller for a problem that in its original formulation is NP- 
hard, as pointed out at the end of Section II. 

The rest of the paper is organized as follows. In Section 
II we describe the system model. In Section III we show 
how message passing can be tailored to solve the problem 
of resource allocation in the uplink of a cellular system. In 
Section IV we present simulation results. Finally, in Section 
V we discuss future work and draw our conclusions. 

II. System model 

We focus on the problem of channel allocation for the uplink 
of an OFDMA system. The overall frequency bandwidth is 
divided into orthogonal sub-carriers and, to reduce allocation 
complexity, we group sets of adjacent subcarriers into F sub- 
channels. As long as the bandwidth spanned by a subchannel 
is smaller than the channel coherence bandwidth, the channel 
spectrum can be approximated as flat in the subchannel. 
Thus, we can assume that the choice of performing resource 
allocation on subchannels rather than on subcarriers causes 
almost no loss in diversity. 

Allocation is performed with the goal of minimizing the 
overall transmitted power subject to rate constraints per user. 
Due to practical considerations, we consider only a limited 
set Q = {0, . . . , Q} of possible transmission formats. A given 
transmission format q corresponds to the usage of a certain 
error correction code and symbol modulation that leads to 
a spectral efficiency r\ q : a user employing format q on a 
certain subchannel transmits with rate R = Brj q , B being 
the bandwidth of each subchannel. The spectral efficiency 

'the algorithm is also known as the belief propagation algorithm 



associated with format q = is n = 0, i.e. no transmission 
at all. The target SNR to achieve the spectral efficiency 
r] q = log 2 (1 + SNR(q)) is SNR(q) = 2"* - 1. Let T be 
the set containing the F available subchannnels. Given the 
format q, the power P n j(q) necessary to user n to transmit 
on subchannel / is computed as 



P nJ (q) = SNR(q) 



BNp 



(1) 



where H n j is the channel gain between user n and the BS on 
the /-th link and N is the power spectral density of the zero- 
mean thermal noise. Channel assignment is exclusive: each 
subchannel can be assigned to only one user and with a just 
a single format. 

Our resource allocation problem is a constrained minimiza- 
tion problem in the vector x = [xi t i, . . . ,xn,f], where the 
variable x n j G Q indicates the modulation format of user n 
on subchannel /. The allocation problem has the following 
general form 

minimize /o(x) (2) 
subject to 

d/(x)<l JeT (Cl) 
Mx)>6„ n = l,...,N (C2) 

ffnW < P m ax,n U=1,...,N (C3) 



Here the objective function /o : T> — ► 
of overall power of the allocation x: 



1+ is the cost in terms 



/o(x) 



N 



E E PnjiXnj) 
n=lfer 



(3) 



the domain V = Q NF is the set of all possible transmission 
formats on all subchannels for all users. The inequality con- 
straints functions df : V — > R + represent the condition of 
exclusive allocation for all subchannels 



d/(x) 



JV 

E 

n=l 



I (Xn,f) 



(4) 



where I(x n j) is 1 if 1 < x n j < Q and otherwise. The 
constraints functions h n : V — > M + enforce that each user 
transmits at least with rate b n 



M x ) = 2^ BT] X 



(5) 



The constraints functions g n : V — > R + enforce that each user 
does not exceed its maximum transmitting power P ma x,n 



5nW 



E 



Pn,f(Xn,f) 



(6) 



The radio resource allocation problem introduced above can 
be shown to be NP-hard by a straightforward reduction from 
the NP-hard problem Multiprocessor Scheduling [8], even 
when only one single transmission format is considered. 



III. Resource Allocation via Message Passing 

In the following, we formulate the allocation problem in 
such a way that can be solved with a message passing 
technique (MP). The advantage of MP is that the computation 
load is distributed among the various nodes by locally passing 
simple messages among simple processors whose operations 
lead, after some time, to the solution of a global problem. 

First of all, to simplify the allocation task we assume 
that each user selects a subset of all available subchannels. 
Let V n C T be the subset of cardinality P < F of 
subchannels that can be allocated to user n. In other 
terms, we assume that x n j may be different from zero 
only if / e V n - As for the choice of the subchannels 
in V n , we make the natural assumption that they 
represent the P best subchannels for user n, i.e. V n = 
{/ € T : l-Hn,/ 1 is one °f the P largest values for user n}. 
Each user may pre-compute its subset of channels before the 
resource allocation algorithm is initiated 2 . 
For our scope, it is convenient to interpret the resource 
allocation problem as a minimum cost problem, where the 
unfulfillment of constraints in (2) gives an infinite cost. 
Thus, we take care of the constraints Cl by introducing the 
cost function C(f) (f G J 7 ), which is if the exclusive 
requirement on subchannel / is fulfilled and oo otherwise 

r o if E i(x n ,f)<i 

C(f) = < neM(f) (7) 



where Af(f) is the subset of users that might use subcarrier 
/, i.e. Af(f) = {n: f G V n }. The constraints C2 and C3 
are dealt by introducing the set of functions W(n) (n = 
1, . . . , N), defined as: 



W{n) 



E Pnj(x„j) 

feVn 



>b n 



oo 



if E Br, x 
feVn 

E Pn,f(Xn,f)<P, 

otherwise 



(8) 

Despite notation complexity, the meaning of (8) is straightfor- 
ward: W(n) is the power transmitted by user n if power and 
rate constraints for user n are fulfilled, and oo otherwise. 
Given the above, it is straightforward to rewrite the resource 
allocation problem in (2) as: 



N 



arg min ( C (f) + E W ^ 
fe r n=i 



(9) 



Since the goal is to get a distributed solution for the above 
minimization problem, we focus on a single variable, e.g., 
x n j, and rewrite the same problem in a form suited for MP 
implementation as: 



x n j — arg mm 

X„,f 



nun 



N 

C(f) + £ W(n) 

n=l 



(10) 



2 We assume perfect channel state estimation between each user and its 
serving BS. 



where notation min denote the minimum over all variables x 
except x n j. 

A. MP implementation 

The MP algorithm has been broadly used in the last years in 
channel coding applications. In particular, when dealing with 
bitwise MAP channel decoding, MP finds an optimum solution 
for the sum-product problem, provided that the correspondent 
factor graph is a tree [9]. The MP algorithm for the sum- 
product problem derives by the distributive law, i.e., by the 
property = II E- However, since min£ = £ min, 

the same property still holds for min-sum problems, where 
minimization replaces addition in the original formulation 3 and 
addition replaces multiplication. By exploiting such a formal 
equivalence, it is straightforward to adapt the MP algorithm 
to the min-sum problem (10). To elaborate, let associate 
with problem (10) a factor graph, where variables Xi yP are 
circular nodes and functions C(f) and W(n) are square nodes. 
Variable nodes are connected with function nodes by an edge 
if and only if the variable appears in the function, i.e. x n , p is 
connected to the P functions C(£) with I G V n and to W(n). 
The factor graph for (10) is depicted in Fig. 1 where we denote 
by x n .p {p = L • • • , P) the transmit format for user n on the 
p-th ordered element of V n . Following a MP strategy, variable 
and function nodes exchange messages along their connecting 
edges until variable nodes can decide on the value of Xi tP . 

Let now assume that the factor graph is a single tree, 
i.e., a connected graph where there is an unique path 
to connect two nodes. In this case, the implementation 
of the MP approach is straightforward. Let firstly intro- 
duce messages as (Q + 1) -dimensional vectors, denoted 
by m = {m(0), m(l), . . . , m(Q)}. In particular, denote 
by m^P/m^j ^ messages exchanged between the C 
function nodes and the connected variable nodes, and by 
m nj V * > 7 m ij^ messages exchanged between the W function 
nodes and the connected variable nodes. 

As in the classical sum-product scenario, message passing 
starts at the leaf nodes, i.e., those nodes which have only 
one connecting edge. In particular, each variable leaf node 
passes an all-zero message to its adjacent function node, whilst 
each function leaf node passes the value of the function to its 
adjacent node. After initialization at leaf nodes, for every node 
we can compute the outgoing message as soon as all incoming 
messages along all other connected nodes are received. 

As far as variable nodes are of concern, the outgoing 
message sent over an edge is simply evaluated by summing 
all messages received from the other edges. With regard to 
function nodes, let first consider the C(f) nodes and focus on 
generic subchannel I. The square node corresponding to C{£) 
is connected to all the variable nodes x n j with n G Af(£). 
If we consider without loss of generality the message to 
be delivered to xj,e (j G M{t)), the q-th element of the 
return message is the solution of the following minimization 

3 See [7], [9] for a detailed description of MP algorithm for the sum-product 
problem. 
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W(l) W(2) W(N) 



Fig. 1. Factor graph for RRM. For ease of representation, we denote by x n ,p 
(p = 1, . . . , P) the transmit format for user n on the p-th ordered element 
ofV„. 



problem: 

m fl V) (l) = min E ™£ r P( x n,t) 

subject to (11) 

E T{x n j)+X(q)<l 
neAT(i),njtj 

In a similar way, we can rewrite message passing rule 
for W(n) nodes. Let focus on generic user u, the square 
node corresponding to W(u) is connected to the variable 
nodes x u j with / G V u . If we consider without loss of 
generality the message to be delivered to x UM (v G V u ), the 
g-th element of the return message is the solution of the 
following minimization problem: 

mr%)=min E P U J (*.,/) + (*„,/) 
fer u 

subject to 

E B Vxu J + Br, q > b u (12) 

feTuJ?" 

E Pu,f(, x u,f*) "F Pum{q) — Pmax.u 

fev u ,f^ 

When a message has been sent in both directions along 
every edge the algorithm stops. It is worth noting that in 
the considered OFDMA cellular scenario the W(n) function 
node and its connected variable nodes are located at the n- 
th user, while all C(f) function nodes are located at the 
BS. Hence, sending messages from variable nodes to C(f) 
function nodes and vice-versa requires actual transmission on 
the radio channel. Instead, message exchange between variable 
nodes and W(n) function nodes is performed locally at the 
users' terminals, without any transmission. 
The solution of Problem (12) requires by far the largest 
computational effort, since it calls for an exhaustive search 
over all possible combinations of transmission formats. Thus, 
in the following we present a new formulation of (12) to find 
the optimal solution with limited complexity. 

B. A Dynamic programming algorithm 

Given a user u, Problem (12) basically consists in finding a 
set of subchannels / G V u , and for each selected subchannel 



the related transmission format to use by u, so that a given 
function is minimized. Such problem can formulated as an 
Integer Linear Programming (ILP) problem introducing binary 
variables y/^ equal to 1 if the user transmits on the subchannel 
/ with the format h, and otherwise. In a general form, such 
a problem can be rewritten as 



f,h 



B VhVf,h > (3 



(13) 



f,h 

^2 p u,f(h)yf,h < « 

$>/,/, < 1 feV u 

h 

y/,h e {0, 1} 

where the cost c/,/, is the cost for user u of transmitting with 
format h on subchannel / (i.e., c/^ = P u j (ft) + (ft)), 
(3 = b u -Br\ q , a = P ma x,u~Pu,v(q)- As in (12), the first two 
constraints correspond to the requirements on the bit-rate b u 
and on the maximal transmission power P ma x,u> respectively. 
The subsequent P constraints impose that at most one format 
is selected for each subchannel /. Note that, b u is limited 
from above by Bt\qP, and such a value is obtained when 
user u transmits with format Q on all the subchannels in P u . 
Assuming that all admissible formats are multiple integer of 
a given spectral efficiency rj, ie r]h = hrj (ft = 0, . . . , Q), we 
can divide all terms of the first constraint of Problem (13) by 
Bf\ to obtain the equivalent constraint 



h y t\h > 

f,h 



Br) 



(14) 



where -JL is limited from above by QP. Observe that, the 
coefficients h in the left-hand side of constraint (14) are integer 
values. Hence, since the left-hand side of (14) is integer, for 
any choice of variables y/^ e {0,1}, -JL can be rounded 



Bf, 



to Li=^J- Moreover, we may assume that values P u j(h) and 
a are integer (e.g., by multiplying all terms of the second 
constraint of Problem (13) by a suitable large number). 

In the following, we show that Problem (13) can be solved 
by a dynamic programming approach [10]. Let z p {d,k) be 
the optimal solution value of Problem (13) defined on the 
first p subchannels, with a "bit-rate" of |_-J^J = d and a 
restricted maximal transmission power of k. We assume that 
z p (d, k) = +oo if no feasible solution exists. Initially we set 
zo(0, k) = and z (d, k) = +oo for all d = 1, . . . , QP and 
k = 0, 1, . . . , a. To compute z p (d, k), we can use the recursion 
(15), where we assume that the minimum operator returns +oo 
if we are minimizing over an empty set. An optimal solution 
of Problem (13) can be found computing zp(QP,a), and 
choosing the minimum of zp(j, a) for j = [^\, ■ ■ ■ ,QP- 
The formula (15) requires the comparison of Q terms, and, 
hence, the optimal solution value of Problem (13) can be 
found in 0(P 2 Q 2 a) operations, only pseudopolynomial, since 



a depends on the input data (i.e, P ma x,u and P u ,u(q))- 
Observe that, if no requirement is given on the maximum 
transmission power used by each user, i.e., if the constraint 
on a can be relaxed, Problem (13) can be solved in 0(P 2 Q 2 ) 
operations, polynomial in the number of subchannels in V u 
and transmission formats. 

C. MP scheduling and Peeling procedure 

As in traditional MP approach for the sum-product prob- 
lem, if the factor graph is a tree there is a natural sched- 
ule for MP given by starting at the leaf nodes and send- 
ing a message once all incoming messages required for 
the computation have arrived [11]. Unfortunately, in general 
the graph which represents minimization problem (10) is 
not a tree (e.g., in Fig. 1 we have a cycle given by the 
path xi.i, C(f),XN,i,XN,p, C(F),xi i p,x\,i). In this case, to 
completely define the algorithm for a generic factor graph we 
need to specify a schedule. It is worth noting that, even if mes- 
sage passing in the presence of cycles is strictly suboptimal, 
the solution found by means of iterative approaches is in most 
cases very close to the optimum (e.g., in the case of bitwise 
MAP decoding of linear block codes) [12], [13]. 
Iterative MP starts at variable nodes, which send an all zero 
message m^ 14 ^ = to their adjacent W(n) function nodes 
and then the algorithm proceeds in iterations. The pseudocode 
of Algorithm 1 illustrates the iterative MP algorithm for a 
generic user n. After I iterations, each user peels off all 
variable nodes > 0. All these variable nodes, say it 

fulfilled nodes, send a message to the BS to communicate that 
the corresponding subchannels have been reserved and the BS 
signals it to all users via the downlink broadcast channel. 

Algorithm 1 Iterative MP procedure for user n 

while J2 Br/ Xn f < b n do 
fev n 

(VW) 



(/ € Vn) 



send m£7° to W(n) (/ € V n ) 

for iter = to I do 

evaluates m { ™ f v) according to (12) (/ € V n ) 



JVC) 

V/ 



send m^p to C(f) (/ £ V n ) 

while not received all m^J^ 1 from C(f) (/ £ V n ) do 

wait 
end while 

end for 

(CV) . (WV) r p — \ 

x nJ <- arg min m nJ (q) (/ G V„) 

9=0,1, ...,Q 

Sends a message containing assignments x n j to the BS 
end while 



(wv) 
V/ 



(/ G Pn) 



Variables corresponding to fulfilled nodes are fixed and do 
not pass any message anymore. At this stage, all users evaluate 
wether they fulfill their rate constraints or not. Those users 
that satisfy their constraints stop participating to MP. All other 
users take part to successive iterations of MP, after having 



z p (d, k) = min < 



z p _i{d, k) + Cpfi (subchannel p is not used by the user) 

— 1, k — P U;P (1)) + c p ,i if d — 1 > and k — P U;P (1) > (p is used with format 1) 



(15) 



z p -i(d — Q,k — P U , P {Q)) + c Pt Q if d — Q > and k — P U , P (Q) > (p is used with format Q) 



updated their rate constraints in (12) on the base of the amount 
of resources they have been allocated. Before starting a new 
cycle of I iterations, each user computes again the set V of 
the best P subchannels among all subchannels which have not 
been yet assigned to other users. The process continues until 
the rate constraint is fulfilled. 

IV. Numerical results 

In this section we present the numerical results of the 
proposed algorithm. We have considered an hexagonal cell 
of radius R = 500 m. The uplink bandwidth is TK = 5 MHz 
so the sampling time is T c = 200 ns. Channel attenuation 
is due to path loss, proportional to the distance between 
the BS and the MS, and fading. The path loss exponent is 
a = 4. We consider a population of data users with very 
limited mobility so that the channel coherence time can be 
assumed very long. The propagation channel is frequency- 
selective Rayleigh fading. The power of the j-th path is: <r| = 

o\ exp ^— , (j = 1, . . . , Nj) where a\ is a normalization 
factor chosen such that the average power of the channel is 
normalized to the value of the path loss, a n = <r T /T c is the 
normalized delay spread with a T = 0.5 [is and Nj = [3a n \ 
is the number of paths taken into an account. 

The available bandwidth W is divided in F — 32 subchan- 
nels and there are N active users at one time. We assume that 
all users request the same rate i.e., b n = b , (n = 1, . . . , N), 
so that bo = Wr] aV g/N, where 7] avg is the average spectral 
efficiency in the cell. The results shown in the following have 
been obtained by setting rj avg = 1 b/s/Hz and averaging on 
500 channel realizations. We compare the performance of the 
proposed MP algorithm with two other resource allocation 
strategies: 

1) The heuristic algorithm presented in [3] that we have 
indicated with the acronym BRCG (Babs + RCG) that 
solves the problem (2) by dividing it in three subprob- 
lems: 1) Decide the number of subcarriers each user 
gets based on rate requirements and the users average 
channel gain (bandwidth assignment based on SNR, 
BABS); 2) Select which subcarriers to allocate to each 
user according a greedy strategy (rate craving greedy, 
RCG); 3) Set the modulation for each subcarriers by 
employing a single-user bit loading technique. 

2) A linear programming (LP) implementation of the allo- 
cation problem (2) formulated as in [4] where the rate 
constraints are translated into a number of subchannels 
to assign to each user. In our implementation, we set a 
unique transmission format for all users on all subchan- 
nels, so that each user is assigned the same number of 
subchannels F/N and transmits with spectral efficiency 



tj = rj avg . By doing so, we neglect on purpose the impact 
of frequency diversity to focus only on the impact of 
multi-user diversity on the allocation performance. 
As far as MP and BRCG are of concern, we set Q = 4, i.e., 
we consider four different transmission formats. Fig. 2 shows 
the average total transmitting power for different number of 
users. As far as MP parameters are of concern, we set P = 
4, 8, 12, 24 for N = 2,4, 8, 16, respectively. Moreover, since 
both BRCG and LP do not take into account any constraint 
on transmitting power, we set maximum transmitting power to 
+oo in (12). Note that the proposed MP algorithm requires the 
minimum average power in all cases. In particular, for small 
values of N, MP and BRCG clearly outperform LP, whilst 
for high values of N, MP and LP outperform BRCG. With 
few users, i.e., for small values of N, each user is assigned a 
great number of subchannels. In this case, the use of multiple 
transmission formats in MP and BRCG algorithms allows the 
transmitter to concentrate the power on the best channels while 
turning off the worst ones. Although designed to achieve an 
optimal global subchannel allocation, the LP scheme shows 
poor performance since it is forced to use the same transmis- 
sion format over all assigned subchannels. On the other hand, 
increasing the number of users reduces the number of channels 
allocated and each user is assigned only 'good' channels, 
thus exploiting the so-called multi-user diversity. Fig. 2 shows 
that for all three algorithms increasing the number of users 
determines a reduction of the average transmission power. 
However, the BRCG algorithm can not fully exploit multi-user 
diversity, since greedy channel assignments are sub-optimal 
and, even using several transmission formats, is outperformed 
by the LP scheme already for N > 8. 

Similar considerations can be drawn when considering the 
outage probability curves, i.e., Figs. 3-4. The outage prob- 
ability P a is evaluated by specifying a maximum allowable 
transmitting power P max for each user. Such a maximum 
power is included in (12), so that outage events in the MP case 
occur when the iterative MP algorithm is not able to provide 
a feasible solution for all users. Differently, since in the LP 
and BRCG approaches we have not included constraints on 
the maximum transmitting power, outage events occur when, 
after the allocation, the power transmitted by a user exceeds 
Pmax- Figs. 3-4 show P as a function of P max for N = 2 
and N = 16. The proposed MP scheme achieves the lower 
outage probability in all cases, thus confirming that it allows 
to perform an optimal subchannel assignment and to profitably 
exploit both frequency and multi-user diversity. Furthermore, 
although the computational complexity of MP depends on 
the number of iterations, at each iteration, the computation 
is naturally distributed among the transmitters, which have to 



2.6 
2.4 

2.2 
2 

S 1-8 
o 

°- 1.6 
d) I 

< 1.4 
12 
1 

0.8 



MP 
-O- BRCG 
-*- LP 



2 4 6 



10 12 14 16 



-12 -10 
P ™, ( dB > 



!■ MP 
o Greedy 

- * - LP 



Fig. 2. Average power consumption versus number of users. 



Fig. 4. Outage probability versus maximum transmitting power for N = 16. 
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Fig. 3. Outage probability versus maximum transmitting power for N = 2. 



solve low-complexity and in practice small problems. 

V. Conclusion 

We have proposed a novel distributed resource allocation 
scheme for the up-link of a cellular multi-carrier system based 
on the message passing (MP) algorithm. Resource allocation 
may rely on a similar MP procedure: with the goal of achieving 
a global optimal assignment, each transmitter iteratively sends 
and receives information messages to/from the base station 
until an allocation decision is taken. The exchanged messages 
are the solution of small distributed allocation problems. To 
reduce the computational load, the MP problems at the termi- 
nals follow a dynamic programming formulation. Hence, even 
if the computational complexity of MP depends on the number 
of iterations, at each iteration, the computation is naturally 
distributed among the transmitters, which have to solve low- 
complexity and in practice small problems. Numerical results 
show that the proposed approach is an excellent solution to 
the resource allocation problem for a single-cell multi-carrier 
system. Moreover, the distributed nature of the proposed 
strategy make it naturally suitable for larger scale resource 



allocation problems, such as global resource optimization in 
multi-cell OFMA systems. 
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